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Abstract 

This contribution demonstrates the potential gain for the quality of results in a simulation of pedes- 
trians when estimated remaining travel time is considered as a determining factor for the movement of 
simulated pedestrians. This is done twice: once for a force-based model and once for a cellular automata- 
based model. The results show that for the (degree of realism of) simulation results it is more relevant 
if estimated remaining travel time is considered or not than which modeling technique is chosen - here 
force-based vs. cellular automata - which normally is considered to be the most basic choice of modeling 
approach. 

1 Introduction 

Users of crowd simulation software tools have noted that these tools mostly produce a too high pedestrian 
density close to a corner when a large group of pedestrians has to walk around this corner compared to real 
people in the same situation [1]. This leads to considerably increased travel times of simulated pedestrians 
compared to real people as well as compared to travel times of simulated pedestrians on a straight path of 
the same length. The same authors express their conviction that real people on contrast organize smarter 
in such a situation. Real people do not insist to walk a short path, but rather accept a longer walking 
path in exchange for a moderate local density and therefore can maintain their walking speed and in final 
consequence achieve a smaller travel time, respectively a travel time which is only slightly increased compared 
to them walking a straight path with the same length. Concerning a first, qualitative empirical evidence, 
someone reading this article in an office in a busy CBD area can witness this behavior by looking out the 
window. For everyone with a less busy window view, YouTube can help [2]. 

Even if one had no empirical indications of human walking behavior in such a situation, it would be a 
good guess that people desire mostly to minimize their travel time and not their travel distance. After all 
traffic planning for road networks rests on the assumption that one of the most important factors - if not the 
most important one - for route choice of vehicle drivers is travel time and not route length. The same people 
who steer vehicles make up for a large share of pedestrians in most situations where pedestrians occur. They 
may move in different circumstances with different implications of their route and speed choices for their 
personal effort, but they are the same people. 

2 Integrating Travel Time in(to) Models of Pedestrian Dynamics 

As real people, when we are in a hurry and desire most to minimize time to arrival, we are faced with 
the problem that we do not know which direction will be the one that allows earliest arrival - imagine for 
example to walk through a busy train station while being late in time for the train. Even in retrospective 
when having arrived it is not always possible to tell if the arrival actually was the earliest possible or if 
along the path one of the many decisions for a movement direction was wrong with regard to minimizing 
travel time. This is not only a consequence of the complexity of the problem, but also simply of incomplete 
information. 



1 



Simulating pedestrians we are in a similar situation. If with given constraints for every point in time, 
every position, and every destination the exact remaining travel time was known for every type of pedestrian, 
it would in fact not be necessary to simulate anymore. For what would be the additional information one 
could gain from such a simulation? The simulation would be more or less a discrete representation or 
instance of the multi-dimensional field of remaining travel times. Furthermore such a simulation would show 
an optimal solution which most probably would not be a realistic one. 

Therefore what is possible is to estimate a map of remaining travel times in each time step and make 
pedestrians choose their desired or preferred walking direction according to this map respectively field. The 
estimation process at first should model realistically the travel time estimation real people do about the 
direction options they have and second it should be self-consistent with the actual model of pedestrian 
dynamics in which the travel time estimation module is used. 

This is not a simple task and as a consequence models of pedestrian dynamics usually set the direction 
of the shortest path as the main or preferred walking direction [3, 4]. However, in the last few years at least 
two such combined systems of a remaining travel time estimation and its usage in a pedestrian simulation 
have been introduced. Both have in common that the travel time estimation (or delay time estimation) is 
first done for small areas and then numerically integrated to receive the field of remaining travel time from a 
spot to the destination, which then is used in the pedestrian model where it amends or replaces (the impact 
of) the field of remaining distance to destination (aka "static potential"). 

The first model is the F.A.S.T. model [5-8]. The F.A.S.T. model has developed out of a cellular automata 
approach. It is arguable if it still is a cellular automata [9], however, the important property that the 
pedestrians move like chess figures on a regular grid is preserved, making the model both: very fast for 
execution [10] and rather coarse-grained in the results. A second property is that it is attractiveness of spots 
which are used in the computation for the next movement step. This distinguishes it from the Social Force 
Model, where the property which is directly calculated is acceleration (and locations follow via integration 
over time). The method to estimate travel times considers only the spatial distribution of pedestrians and 
not their dynamic movement properties 1 [11]. The method for numerical integration is the so called variant 
2 as introduced in [12]. Compared to the Euclidean metric - if a field of constant local distances or travel 
times is integrated - it includes some artifacts, but the way it is used in the pedestrian model reduces them 
largely and the computation time of variant 2 is below the exact methods mentioned in the next paragraph. 
In this way a main intention of the F.A.S.T. model - having a very fast pedestrian simulation - is preserved, 
even when the travel time estimation module is active. 

The second model is the Social Force Model [13, 14] respectively its implementation in PTV Viswalk 
[15, 16]. In this approach of a combination of a pedestrian dynamics model with a travel time estimation 
the velocity of a pedestrian is considered when it is estimated which delay time he will impose to another 
pedestrian in his surrounding [17] 2 . For the numerical integration the Fast Iterative or Fast Marching Method 
is used. Either method - in case of distance computation - yields the nearly correct Euclidean distances 
[18-21]. In the Social Force Model the negative gradients of the resulting field of estimated remaining travel 
time are used as direction of the desired velocity of a pedestrian being located at the field's corresponding 
position. Apart from that PTV Viswalk's Social Force Model remains unchanged. 

Both models have been used with their travel time estimation module in different use cases [22-27]. 
However, in this contribution it is not intended to solely demonstrate the effect of the usage of the field 
of remaining travel times, but to argue that the decision to use it or not is more important (i.e. has a 
greater effect) than the decision which basic modelling approach - force-based or cellular automata-based - 
is applied. 

2.1 Implications and Interpretation 

Normally this would be the place to define the math of the travel time estimation module(s). However, this 
has been done extensively before [11, 17]. Instead of a repetition this subsection discusses some implications 
of the method, i.e. it takes some interpretative steps. 

1 Of course obstacles are also considered. Considering only them would be no big deal, as they are static. The field of 
remaining travel times would then be identical or at least very similar to the field of remaining distances. And it would only 
be necessary to compute it once in advance of the simulation. One of the difficulties of considering travel time is that the field 
has to be recalculated in every time step or at least frequently. 

2 Also see that paper for an extensive discussion of related work. 
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It has been observed [28, 29] and sometimes criticized [30] that the Social Force Model (and most if not 
all of its variants) strictly superpose the repulsive forces of pedestrians. Superposition means that the effect 
of any pedestrian A on pedestrian B is absolutely independent of any other pedestrian in a scenario. 

Superposition is easiest noticed for the Social Force Model due to the similarities with Newtonian Me- 
chanics, however, if one thinks of superposition of effects not only forces, then many models with a different 
modeling approach are just as well superposing effects of pedestrians. 

In physics superposition is so common that one might even wonder, how else a model should be con- 
structed. One approach to the answer is to formulate the extreme opposite: "the effect of any pedestrian 
Ai onto a specific pedestrian B is depends on all other pedestrians Aj and their states". This phrase is a 
concretion for (the simulation of) pedestrians of Sherif's general statement on social systems: "the properties 
of any part are determined by its membership in the total functional system" [31] and it is a holistic position. 
Figure 1 shows a simple situation, where a superposition of effects does not yield a fully realistic effect and 
where therefore applying a holisitc approach improves the results. 




Figure 1: An inequality with regard to the effects of the displayed situations. The idea is that the observer 
here basically wants to move straight on, but is faced with pedestrians walking the opposite direction. The 
two white pedestrians alone would eventually "channel" the observer to pass between them. In a force-based 
model their deviation triggering effects - at least for this very moment - would largely cancel. The red 
pedestrian alone would make the observer do a small detour to the left or to the right. However, together 
the three trigger the observer to do a larger detour to the left or to the right. The effect of all three - 
especially the red one - on the observer is determined by their membership in the total functional system. 

The module of travel time estimation by construction is a holistic approach. One single pedestrian 
either can have a negligible effect in one setting and a very large one in a different setting, depending on 
the distribution and dynamic properties of all other pedestrians (and obstacles) in the system, even if his 
position and dynamic properties are the same in both settings. Therefore the holistic travel time estimation 
module is an ideal supplement to the superposing force approach. 

This line of argumentation shall not mislead to the conclusion that a holistic (concrete: travel time 
estimating) approach alone is sufficient or superior to a superposing (concrete: force-based) approach. There 
are situations which are better addressed with forces. Imagine for example in Figure 1 that there are walls 
to the left and the right of the group of three pedestrians. Then the observer would be forced to slow down, 
even be forced to stop. This is better modeled with forces. The travel time estimation module then is not of 
much help, as it is there to compute a good walking direction. Confined between two walls there would be 
no more really good direction, but just the task to slow down and resolve the situation without collisions. 

A different example which fits better for force-based modeling is when the density is high enough that 
each pedestrian is forced to go with the flow. 
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On the contrary scenarios which are only sparsely populated with individual pedestrians can be addressed 
with both approaches alike. 

It is possible to some extend to use only one of the two approaches to model effects which are better 
modeled with the other approach (only forces and a simple determination of the desired walking direction 
as has been done usually or only a holistic approach and only a basic treatment of forces and accelerations 
[32]), but it makes life easier if both are applied, each for its own "natural" purpose. 

With the emphasis in the previous lines it has already been insinuated: in a nutshell it appears that it 
is more natural and therefore easier to use forces to model behavior which is experienced to be forced 3 . A 
main reason for this experience often is the small time available to react on a stimulus which does not allow 
for decent planning. On the contrary the holistic approach with its potentially unlimited range for strong 
interactions is better suited to model conscious planning, intentions or desires. 

3 The Effect of Integrating Travel Time into Models of Pedestrian 
Dynamics 

3.1 Example 1 

The scenario as shown in Figure 2 is simulated both with the F.A.S.T. model and the Social Force Model of 
PTV Viswalk without applying their modules of estimated remaining travel time. The geometry mimics a 
simplified part of a stadium and its exterior. In the beginning two groups with each 5,000 pedestrians start 
their movement in the interior of the stadium. 




Figure 2: Scenario. Pedestrians can move on light grey areas, while dark grey areas denote obstacles. 
Pedestrians are set on the light red and light blue areas into the scenario and have to move to the green area 
at the upper right. 

The first observation is that it is possible to find parameters for both models with which the simulation 
results become very similar first for the total evacuation times, second for the average individual egress 
time and third for the spatial en gros distribution of pedestrians all throughout the simulation. If one 
has a detailed look at the trajectories, these look differently for both models. However, the idea here is 
demonstrating equivalence for someone interested in an overview, i.e. envisioned as "looking from high 
above on the whole scenario" . 

3 Note that not in every language the same word is used for physical force and an (en)forced action. In German it is for 
example "Kraft" and "Zwang". 
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These are the parameters that were used with the F.A.S.T. model: ks — 1.0, hp = 0.0, ki = kw — 0.6, 
\i = 0.0; and these with PTV Viswalk: r = 0.41, A socMo = 1.0, B socMo = 0.2, A = 0.0, A soc ^ mean = 0.4, 
Bsoc,mean = 2.8, VD = 1.7, noise = 1.2. The latter parameters are identical or close to the default values 
except for A soc ^ so , whose value has been clearly reduced to achieve the relatively high densities which occur 
in the F.A.S.T model with a comparable regularity. 

In both cases the free (or desired or preferred) speeds have been (nearly) equally distributed between 1.2 
and 2.0 m/s. As they are defined with different methods in both models, the distribution of free speeds is 
be not exactly identical, however very similar. Note that it has been found to be crucial for the equivalence 
of the results of both approaches to have a very similar distribution of free speeds. 

Someone familiar with one of the models will immediately see that neither the chosen parameters nor 
the free speeds are extreme in any way. 

With these settings the individual egress and the total evacuation times result as shown in Table 1. 





Av. Individual Egress Time 
Min., Average ± STD, Max. 


Total Evacuation Time 
Min., Average ± STD, Max. 


F.A.S.T. 
Viswalk 


1317.1 1332.2 ± 4.9 1347.9 
1324.0 1337.4 ± 4.6 1350.6 


2529.0 2558.9 ± 9.3 2590.0 
2528.2 2555.9 ± 9.1 2584.3 



Table 1: Results statistics from 1,000 simulation runs with each of the models. "Av. Individual Egress 
Time" means the average of egress times of all pedestrians of one simulation run. Over 1,000 runs there are 
1,000 averages of which the minimum, average, and maximum are displayed. 

Selecting for both models the simulation run which has the least sum of mean square deviations from 
the averages for total evacuation and individual egress times Figure 3 in this way compares the average 
simulation runs of both models and shows that the results are very similar all through time and space. 

While these snapshots show how similar the results are, a second look reveals that the pedestrians move 
unrealistically close around the corners. As a consequence they ignore all but two or three of the exit gates 
shortly before the destination area. One may therefore expect that the times as shown in Table 1 are much 
larger than would be realistic. Someone familiar with one or both of the models (better: model classes) will 
know that this cannot be fixed by just modifying the parameters. Of course simulated egress and evacuation 
times can be reduced to some degree with different parameters, but only at the cost of having the models 
yield unrealistic results somewhere else in this geometry and even more so in different geometries. A reader 
not familiar with these models here needs to trust the author that this is in fact the case as demonstrating 
this would by far exceed the available space of this contribution and a strict formal prove is simply not 
possible. 

Using the module of remaining travel time estimation on contrast greatly reduces egress and evacuation 
times and strongly influences the walking behavior in both models. For these simulations the same parameters 
as above have been used and in addition for the travel time estimation module used with the F.A.S.T. model 
it has been set ksdyn = 1-0 and s a dd = 10 (see [11, 22]) and for the one of PTV Viswalk g = 1.5, h = 0.7 
and the impact strength has been set to 100% (see [17]). For this second stage of simulations it has not been 
attempted actively to find parameters with which the results of the two approaches match, i.e. these values 
are something like the defaults of the methods judged by previous experience. 





Av. Individual Egress Time 
Min., Average ± STD, Max. 


Total Evacuation Time 
Min., Average ± STD, Max. 


F.A.S.T. 
Viswalk 


784.6 787.6 ± 1.8 789.9 
666.5 672.0 ± 2.4 676.2 


1345.0 1352.7 ± 4.9 1363.0 
1227.4 1234.9 ± 4.6 1244.2 



Table 2: Results statistics from 15 simulation runs with each of the models when the travel time estimation 
module is used. 

Figure 4 shows snapshots from both simulations. By comparing with Figure 3 it can be seen that the 
travel time estimation module modifies the simulation dynamics of both models into the same direction: 
pedestrians utilize the available space at corners better, they utilize the full width of the long corridor to the 
right better and they utilize the exit gates better compared to the simulations without travel time estimation 
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method as shown in Figure 3 and the red pedestrians - respectively a center line through the stream of red 
pedestrians - in the open area do not head for the corner of the obstacle but (about) for the edge of the 
group of blue pedestrians. Last but not least: for both models the evacuation and egress times are about 
halved as a comparison of Tables 1 and 2 shows. 




Figure 4: Snapshots from the simulations with dynamic potential: left: F.A.S.T., right: Viswalk, each at 
time 350 seconds. 

3.2 Example 2 

The first example was composed of a number of basic elements: corridors, corners, gates. To clarify the effect 
of the travel time modules this is now to be amended by a most simple example: a group of pedestrians 
walks around a corner and this is compared with examples where the same number of pedestrians walks 
comparable distances straight on. Details of the geometry are shown in Figure 5 

The corner situation is simulated once without and once with activated travel time estimation module. 
For the straight examples the travel time estimation is not used it would not have much effect there anyway. 
So why simulate it at all? The straight geometries are only there to set a reference frame which makes the 
results less dependent on specific parameter settings. The value of the simulations in the straight geometries 
will become apparent with the results. 

The results for the F.A.S.T. model are given in Table 3 and for PTV Viswalk in Table 4. In this second 
example it was not intended to achieve agreement between both models. Instead PTV Viswalk was run 
with default parameters and for the desired speed a population composition was chosen as in the RiMEA 
[33] test cases (50% male, 50% female; 4% mobility impaired and beyond that the age groups (<30, 30-50, 
and >50) with equal ratios) with their desired walking speeds according to [34]. In this way much smaller 
desired speeds were involved than for the F.A.S.T. model where preferred speeds between 1.2 and 2.0 m/s 
were used. The other parameters were set as in example 1. 





Individual Egress Time 
Min., Average ± STD, Max. 


Total Evacuation Time 
Min., Average ± STD, Max. 


Corner w/o 
Corner w 
Straight short 
Straight long 


79.3 84.0 ± 1.2 88.6 
66.2 69.4 ± 0.8 72.6 
41.8 43.4 ± 0.5 45.0 
67.2 69.8 ± 0.8 72.5 


142.0 154.2 ± 3.4 171.0 
109.0 114.8 ± 2.3 127.0 
65.0 70.1 ± 1.9 80.0 
102.0 107.4 ± 2.2 118.0 



Table 3: Results statistics from 3000 simulation runs with the F.A.S.T. model. The first line gives the times 
when no travel time estimation module is used, the second when it is. 

Although two different models, two different methods for travel time estimation, and different distribu- 
tions of the desired speed were used, and although for this reason the total evacuation time and the individual 
average egress time differ in their absolute value for F.A.S.T. and PTV Viswalk, the results shown in Tables 
3 and 4 are similar with regard to a number of properties (numbers with regard to average values): 
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Figure 5: The corner and the two straight geometries for comparison. 500 pedestrians have to walk from 
the dark red to the dark green area. The red line is 38 m (28 + 10) long, the cyan one 69.4 m (28 + 10 tt + 
10) long. The width of the areas is 20 m. The dark red insertion area has a length of 6 m. 





Individual Egress Time 
Min., Average ± STD, Max. 


Total Evacuation Time 
Min., Averagei STD, Max. 


Corner w/o 
Corner w 
Straight short 
Straight long 


104.2 110.5 ± 1.9 117.0 
78.2 82.5 ± 1.4 88.3 
52.0 55.1 ± 0.9 58.3 
81.9 86.6 ± 1.4 90.9 


229.4 248.1 ± 6.6 274.4 
157.9 174.2 ± 5.0 191.9 
96.6 111.1 ± 3.4 123.2 
150.4 176.2 ± 5.4 188.8 



Table 4: Results statistics from 3000 simulation runs with PTV Viswalk. The first line gives the times when 
no travel time estimation module is used, the second when it is. 
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- The individual average egress time in the corner geometry is 120% (F.A.S.T) respectively 128% (PTV 
Viswalk) compared to the case of the straight long corridor if no travel time estimation is made. 
However, it drops to 99% respectively 95%, if the travel time estimation module is active. Simpler 
spoken: they are clearly longer without and similar with travel time estimation. 

- For the total evacuation time about the same relation holds: there the time drops from 144% (F.A.S.T) 
respectively 141% (PTV Viswalk) to 107% respectively 99%, if the travel time estimation module is 
active. Again the summary is that the time is clearly longer in the corner geometry compared to the 
straight long corridor if the travel time estimation module is not used, but it is about the same with 
both models if it is used. 

- With the corner geometry the average individual egress time is 83% (F.A.S.T.) respectively 75% (PTV 
Viswalk) when the travel time estimation module is used compared to if it is not used. 

- With the corner geometry the average total evacuation time is 74% respectively 70% when the travel 
time estimation module is used compared to if it is not used. 

- If the travel time estimation module is active the individual average egress time in the corner geometry 
is 160% (F.A.S.T.) respectively 150% (PTV Viswalk) of the one in the straight short corridor. 

- If the travel time estimation module is active the total egress time in the corner geometry is 164% 
(F.A.S.T.) respectively 157% (PTV Viswalk) of the one in the straight short corridor. 

As the properties for F.A.S.T. and PTV Viswalk in all cases are separate by a maximum of ten percent 
points, while the effect of the travel time estimation module is much larger, these relations show that the 
details of the pedestrian model, its parameters (especially walking speeds) and the details of the travel time 
estimation module method and its parameters are of minor relevance. The effect is clear and distinguishes 
the case when it is used from the case when it is not used. This implies that with empirical data it should 
be clearly possible to decide if real people move in such a situation rather based on a shortest or a quickest 
path paradigm, i.e. if it enhances the degree of realism when the travel time estimation module is used. 
In addition it should be possible to compute those relations above where times at the corner geometry are 
related with times in one of the straight corridors with any pedestrian simulation model available. If then 
for example total evacuation time at the corner is about 150% of the one in the straight long corridor it 
can be concluded that pedestrians in the model navigate according to a shortest path paradigm, while if the 
times are about equal or the time in the corner geometry is even shorter than the one in the long straight 
corridor it can be conclued that they follow a quickest path paradigm. 

For completeness Figures 6 and 7 show snapshots from the simulations. 




Figure 6: Snapshots from the simulations with the F.A.S.T. model at 70 seconds. 
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Figure 7: Snapshots from the simulations with Viswalk at 70 seconds. 



4 Summary 

In this contribution it was argued that it is of greater relevance for the degree of simulation results of a 
pedestrian simulation model that pedestrians walk into the direction of estimated earliest arrival instead of 
the direction of the shortest path than the details of a) the model of pedestrian dynamics and b) the details 
of the travel time estimation module. For this it was first shown that in a non-trivial example two different 
approaches to modeling pedestrian dynamics can yield very similar results on a moderate level of aggregation 
if one sets the model parameters appropriately. It was then shown that for both approaches the degree of 
realism is increased, if one changes from a shortest path direction to a quickest path direction paradigm. 

The examination of a simple example scenario confirmed this. With two different models and different 
desired walking speed distributions and different travel time estimation modules it was found that the effect 
of using a travel time estimation module (or using it not) had about the same impact on certain relative 
properties. This implies a comfortable situation for comparison with empirical data. 

A first obvious conclusion is that it would be very interesting to have empirical data for situations as 
the ones in example 2: where between the results for the two straight corridors do the results of the corner 
geometry lie for real people? 

The second conclusion is that because it is a) possible to have some cellular automata-based model yield 
such similar results and b) not possible for both modeling approaches to enhance the degree of realism in 
certain situations, variants of the Social Force Model which modify exclusively or mainly the forces and do 
not include a travel time estimation module [29, 35, 36] most probably cannot help to improve the realism 
of results to the same degree as the travel time estimation module can. 

Finally one can summarize that in the formulation of a pedestrian simulation model we can let guide 
ourselves by a quote which is ascribed to Blaise Pascal: "Desire and 4 Force are responsible for all our 
actions. " A general way to interpret it in this context is that one equally has to take care for desired 
(planned or conscious) and for enforced (reactive, unplanned) determinants of motion. A more concrete way 
to interpret it for a force-based approach to pedestrian dynamics is that one has to put equally much effort 
in modeling the desired velocity as in modeling the forces. 

4 Emphasize added by the author of the facing work. 
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